An efficient algorithm to find a set of nearest elements in a mesh 



o 

(N 
o 

o 

(N 



00 

q 

O 



> 

o 

oo 
O 

d 



X 



Gleb Novitchkov 

email: gleb. novitchkovOgmail. com 

October 4, 2011 



Contents 

1 Introduction 

2 Definition 

2.1 Mesh and the representation of the elements 

3 The algorithm 

4 Parallelization 

5 Conclusion 



1 Introduction 

Here we present an algorithm that find a list of elements neighboring some given element in a 
linear time. More precisely, if there are A/'eiem elements in the mesh, the runtime of the algorithm 

is 0(A^elem). 

2 Definition 

Definition 2.1. By element we mean a 3-simplex imbedded in M^. 
Essentially an element is a tetrahedron in M^. 

Goal: given an element of the mesh, we want to find a set of elements that are near this 
element. 

Definition 2.2. Given an element E, we call an element F a vertex-near element of E if E 
and F share common vertex; we call an element F a edge-near element of E if E and F share 
common edge; we call an element F a face-near element of -E if and F share common face. 
An element F is near element oi E ii F is either face-near, edge-near, of vertex-near element of 
E. 



2.1 Mesh and the representation of the elements 

A mesh is specified by the cloud of n points, or nodes, {po,Pi, ■ ■ ■ ,Pn-i}, where pi = {xi,yi, Zi 
An element Ei is specified by the four nodes, Ei = {piQ,Pij^,Pi2,Pi3}- 
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3 The algorithm 

The idea of the algorithm is to do the histogramming of the nodes. 

Step 1. (initialization) For A^node nodes allocate array L of lists, the length of the array is A^node- 

Step 2. (histogramming) For each element ieiem, < ieiem < -^eiem — 1, do: 

Step 2.1 For each node Pj^, < i < 3 of element icicm ^-dd index icicm to the lists 
L[jo], L[ji], L[j2], L[js]. In C++ notation, it should be: L [jq] .push_back(ieiem) , 
L[ji] .push_back(icicm), L[j2] .push_back(ieiem) , LLjs] .push_back(icicm) • 

Step 3. (Finding neighboring elements) Given element Ei = {jo, ji, 32, js), list of elements neigh- 
boring Ei is given by the union of the lists L[jo], L[ji], L[j2], and L[js]. 

Step 4. Deallocate array L of lists. 

4 Parallelization 

Algorithm admits easy parallelization by OpenMP, one only has to take care to use critical section 
for Step 2.1: adding element index to the lists should be done inside the critical section. 

5 Conclusion 



Missing details will be provided later. 
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